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Statistische Testverf ahren ziim objektiven Nachweis von 
Auditory Steady- State Responses (ASSR) im Frequenzbereich 

Die Erfindung betrifft das Gebiet der objektiven Messung des 
5 Horvermogens unter Nutzung von auditorisch evozierten steady- 
state responses (auditory steady-state responses, ASSR) . Die 
vorgeschlagenen Verfahren zum objektiven Nachweis von ASSR im 
Frequenzbereich konnen sowohl bei einem ASSR-basierten 
Neugeborenen-Horscreening als auch bei der objektiven 
10 Horschwellenmessung auf ASSR-Basis eingesetzt werden* 

Es sind zwei verschiedene Arten von ASSR bekannt : 

1. click-evozierte ASSR 

2. ASSR, die durch einen amplituden- oder f requenzmodulierten 
15 Dauerton evoziert werden, die auch als amplitude-modulation 

following response (AMFR) bezeichnet werden. 
Beide Arten von ASSR werden im Frequenzbereich durch sog. 
Harmonische (eine Grundwelle und mehrere Oberwellen) 
beschrieben (Abbildung 1) . Die Frequenz der Grundwelle 

20 entspricht der Click-Reizrate bzw. der Frequenz des 

Modulationssignals . Die Frequenzen der Oberwellen ergeben 
sich als Vielfache der Frequenz der Grundwelle. Das bedeutet, 
da(i die gesamte Antwort durch wenige Spektrallinien 
reprasentiert wird. Der wesentliche Anteil der durch das 

2 5 Spont an- Elekt roenzephalogramm ( Spent an- EEG) bedingten 

Rauschleistung ist dagegen auf die zwischen den Harmonischen 
liegenden Spektrallinien konzentriert . 

Der objektive Nachweis der bekannten ASSR erfolgt nahezu 
30 ausschlielilich im Frequenzbereich. Im Frequenzbereich ist 

jede Spektrallinie durch eine spektrale Amplitude und einen 
Phasenwinkel definiert. Fiir den Antwort-Nachweis im 
Spektralbereich sind mehrere statistische Verfahren bekannt. 
(Stapells DR, Makeig S, Galambos R. Auditory steady-state 
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responses: Threshold prediction using phase coherence. 
Electroencephalography and Clinical Neurophysiology 
1987;67:260-270; Valdes JL, Perez-Abalo MC, Martin V, Savio 
G, Sierra C, Rodriguez E, Lins O. Comparison of statistical 
5 indicators for the automatic detection of 8 0 Hz auditory 
steady state response (AMFR) . Ear and Hearing 1997;18:420- 
429), die als sog. One-Sample Tests nur den Phasenwinkel . oder 
auch Phase und Amplitude einer einzelnen Spektrallinie, 
vorzugsweise der 1. Harmonischen (Grundwelle) auswerten. 
10 Hierzu ward das registrierte Zeitsignal epochenweise in den 
Frequenzbereich transf ormiert . Die Lange der transf ormierten 
Epochen mu/5 so gewahlt werden, dal3 die Epoche exakt ein 
ganzzahliges Vielfaches der Periodendauer der Click-Reizrate 
bzw. der Modulationsf requenz ist. In dem nach der 
15 Transformation vorliegenden Frequenzspektrum wird die der 
Click-Reizrate bzw, der Modulationsf requenz entsprechende 
Spektrallinie (Grundf requenz) aufgesucht und getestet. Der 
Vorteil des Antwort-Nachweises im Frequenzbereich gegeniiber 
dem direkten Nachweis im Zeitbereich besteht darin, dali im 
20 Spektralbereich der Anteil der Rauschleistung, der durch die 
zwischen den Harmonischen liegenden Spektrallinien 
reprasentiert wird, die Response Detection nicht stort, da 
diese Spektrallinien nicht in die Testung einbezogen werden. 
Der Nachteil der in den oben genannten Publikat ionen 
25 eingesetzten sogenannten "one-sample Tests" besteht in der 

Beschrankung der statistischen Testung auf die Grundf requenz . 

werden jedoch in der Regel nicht allein durch die der 
Click-Folgef requenz bzw. der Modulationsf requenz 
entsprechende Grundf requenz , sondern auch durch eine oder 
30 mehrere Oberwellen reprasentiert, auf die ein nicht zu 
vernachlassigender Anteil der Response-Signalleistung 
entfallt. Ein objektives Nachweisverf ahren, das sich nur auf 
die Grundwelle beschrankt, ist deshalb nicht optimal. 
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Von StCirzebecher sowie von Stiirzebecher et al. werden 
ebenfalls im Frequenzbereich arbeitende statistische Test- 
Verfahren beschrieben (Sturzebecher E: "Method for hearing 
screening of newborn by means of steady- state response evoked 
5 with high click rate", European patent application EP 

01610060,4.; Sturzebecher, Cebulla M, Baag M, Thie R: 
"Verfahren zur objektiven f requenzspezif ischen 
Horschwellenbestimmung mittels der Amplitude-Modulation 
Following Response (AMFR) " , European patent application 

10 EP1099408 A2) . Als statistische Tests werden sogenannte 

q-sample Tests eingesetzt/ die fiir den statistischen Nachweis 
der ASSR sowohl die Grundwelle als auch die relevanten 
Oberwellen einbeziehen. Es handelt sich dabei um den aus der 
Literatur bekannten "q-sample uniform scores Test" (Mardia . 

15 KV, Statistics of directional data. Academic Press London and 
New York 1972) und eine von Sturzebecher et al. 
vorgeschlagene Modifikation dieses Tests (Starzebecher, E^ 
Cebulla M, Baag M, Thie R: "Verfahren zur objektiven 
f requenzspezif ischen Horschwellenbestimmung mittels der 

20 Amplitude-Modulation Following Response (AMFR) " , European 
patent application EP1099408 A2) . 

Die Rechenvorschrif t fur den von Mardia, 1972 beschriebenen 
q-sample uniform scores Test lautet : 

25 Let {xik; l<i<m, l<k<q}bea collection of random 
variables (phase angles cp^-^ ) ; q is the number of samples 
(spectral lines) with the sample size m (number of epochs)., 
i.e. there are q x m = n phase angle values. 

The n phase values were ranked in a single seq[uence . Let rik, 
30 i = 1,..., m, be the ranks of the phase angles in the kth 
sample . 

The phase angles cp^ are then replaced by the uniform scores 
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10 


2-71 


n 

The test statistics used is 


with Q = J]cosP,jt; -S*^ =Ssinp,-^ 


j=i i=i 


where 


rik are the ranks of the n phase angles (n = q x m) , 
q is the number of samples (number of included spectral 
15 lines) 

and m is the sample size (number of epochs) . 

W is distributed as Chi-square with 2 (q-1) degrees of 

freedom. 


20 Wie aus der Rechenvorschrif t erkennbar, werden nur die 
Phasenwinkel genutzt, die spektralen Amplituden bleiben 
unberucksichtigt - Zu diesem Inf ormationsverlust kommt noch 
ein weiterer hinzu. Es gehen nicht die Phasenwinkel selber, 
sondern nur die Ranks in die Berechnung des Testwertes ein. 

25 Das hat zwar den Vorteil, dali das Verfahren 

verteilungsunabhangig (nonparametric) ist, die Folge der in 
Kauf genommenen Inf ormationsverluste ist aber eine geringere 
Power des Tests, Fur das Horscreening und fur die objektive 
Horschwellenbestimmung sollte jedoch die Testpower so hoch 

30 wie moglich sein. Sturzebecher et al. haben deshalb eine 
Modifikation des Tests vorgenommen (hier als Test- 
Modif ikation 1 bezeichnet) , die neben den Phasen auch die 
spektralen Amplituden in Form der Ranks der Amplituden 
berucksichtigt : 
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Test-Modif ikation 1 

Additionally to the phase angles, the spectral amplitudes Ai^ 
were taken into account. Like the phase angles, the spectral 
5 amplitudes Aik are ranked in a single sequence: Let aik, 

i = 1,..., m be the ranks of the spectral amplitudes Aik in 

the kth sample. The phase angles cp,.^ were replaced by the 

uniform scores yg^.^ _ 2 7i r^^ ^ 

n 

10 The test statistics used for the modified q-sample uniform 
scores test is 


15 


with C; =£;a,.^-cosp,;t; '5*=Z«/Asinp,.4 


and ;ff«=^-^. 

n 

20 

where 

rik are the ranks of the n phase angles (n = q x m) , 
q is the number of samples 
25 m is the sample size of the q samples and 

aik are the ranks of the n spectral amplitudes Aijc. 

Diese bereits bekannte Testmodif ikation nutzt zwar durch die 
3 0 Einbeziehung der spektralen Amplituden mehr Information aus 
als der q-sample uniform scores Test von Mardia (1972) und 
ist wegen der Verwendung der Ranks statt der realen Phasen- 
und Amplitudenwerte weiterhin parameterf rei , die Arbeit mit 
den Ranks der Phasen und Amplituden anstelle der 
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tatsachlichen Werte bedeutet aber weiterhin, daft die 
vorhandene Information nicht vollstandig genutzt wird mit der 
Konsequenz, daU die Testpower nicht optimal ist. 

Aufgabe der Erfindung ist es, verschiedene weitere 
Modifikationen des bekannten q- sample uniform scores Test 
vorzustellen, bei denen die zu Verftigung stehende Information 
der Phasen bzw. der Phasen und der spektralen Amplituden 
vollstandig ausgenutzt wird. 

Folgende neue Modifikationen des bekannten q-sample uniform 
scores Tests werden vorgeschlagen: 

Test-Modif ikation 2 
15 Es werden nur die Phasenwinkeln verwendet, aber im 

Unterschied zum bekannten q-sample uniform scores Test von 
Mardia, 1972 wird nicht mit den Ranks, sondern mit den uber 
die Fouriertransformation berechneten Phasenwinkeln 
gearbeitet . 
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25 


30 


Let {xi^; 1 < i < m, 1-^ k ^ q} be a collection of random 
variables (phase angles 9,,); q is the number of samples 
(spectral lines) with the sample size m (number of epochs), 
i.e. there are q x m = n phase angle values. 

The test statistics used is 


^ *=i 


where 

35 
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q is the number of samples (number of included spectral 
lines) 

and m is the sample size (number of epochs) , 

5 

Test -Modif ikation 3 

Es werden wie bei der bekannten Modif ikation spektrale 
Amplituden und Phasenwinkel verwendet, die Phasenwinkel aber 
ungerankt, wahrend fur die spektralen Amplituden weiterhin 
10 die Ranks in den Test eingehen. 

The spectral amplitudes Ai^ are ranked in a single sequence: 
Let aiicf i = 1,..., m be the ranks of the spectral amplitudes 
Aik in the kth sample. 
The test statistics used is 


15 



m 


m 


with 


20 


where 


25 


q is the number of samples 

m is the sample size of the q samples and 

aik are the ranks of the spectral amplitudes Ai^ 

(p^.^ are the phase angles. 
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Test -Modif ikation 4 


Hier werden sowohl die mittels der Fouriertransf omation 
berechneten VJerte der Phasenwinkel als auch der spektralen 
Amplituden direkt (ungerankt) verwendet. 
The test statistics used is 
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92.(^ + 1)2 m 

5 with Cl=Y,Aik'^os(Pi^; 5;=X!4*-sincp,-^ 

1=1 1=1 

where 

10 q is the number; of sainples 

tn is the sample size of the q samples and 
Aik are the spectral amplitudes 
<p^^ are the phase angles. 

15 

Ein wesentliches Problem bei der Anwendung der 
erf indungsgemalien Test-Modif ikationen besteht darin, daB 
infolge der durchgef iihrten Modif ikationen die jeweils 
zugehorige Dichtef unktion der Testwerte (probability density 
20 function) fiir die Nullhypothese unbekannt ist. Deshalb konnen 
die fur die Anwendung der Tests erf orderlichen kritischen 
Testwerte ( er f orderlich fiir die Entscheidung : Testergebnis 
positiv Oder negativ) nicht den in der Literatur vorliegenden 
Tabellen entnommen werden. 

25 

Eine bekannte Moglichkeit fiir die Berechnung der 
Dichtef unktion der Nullhypothese bietet die Monte Carlo 
Simulation.. Hierbei wird..mit einem. Zuf allsgenerator eine sehr 
grofie Anzahl Paare von Zuf allszahlen erzeugt. Aus jedem 

3 0 Zahlen-Paar wird eine spektrale Amplitude und ein 

Phasenwinkel berechnet und darauf der statistische Test 
angewendet* Aus der resultierenden sehr grofien Anzahl von 
Testwerten wird die Verteilung der Nullhypothese berechnet. 
Aus der Verteilung kann der gesuchte kritische Testwert 

35 abgelesen werden. 


. «»i I* 

Die bei den Simulationen angenommene No rmal vert ei lung der 
spektralen Amplituden und Phasen kann jedoch bei den realen 
spektralen Amplituden und Phasenwinkeln nicht vorausgesetzt 
werden . 

5 

Bestandteil der Erfindung ist deshalb das folgende Verfahren 
zur Berechnung der Verteilung der Nullhypothese, das die 
reale Verteilung der spektralen Amplituden und Phasenwinkel 
berucksichtigt : 

10 Voraussetzung ist das Vorliegen einer grofien Anzahl (>100) 
ASSR-Recordings, bei denen die Rohdaten des abgeleiteten 
Elektroenzephalogramms (EEG) kontinuierlich auf der 
Festplatte gespeichert wurden (ca. 200 Epochen^ Lange einer 
Epoche ca. 1 Sekunde) . Wie Abbildung 1 zeigt, ist die 

15 Antwort auf wenige Spektrallinien (Grundwelle bei 160 Hz, die 
hoheren Harmonischen bei Vielfachen von 160 Hz) begrenzt . Da 
bei einer EpochenlSnge von etwa 1 Sekunde die spektrale 
Auflosung etwa 1 Hz betragt^ liegen zwischen zwei 
Harmonischen mehr als 150 Spektrallinien, die nur das aus dem 

20 Spontan-EEG resultierende Rauschen enthalten. Wenn man den.' 
statistischen Test auf diese Spektrallinien anwendet, so 
erhalt man bei 100 Recording von je 200 Epochen etwa 3 000 
000 Testwerte (100x200x150) . Die aus diesen Testwerten 
berechnete Verteilung stellt eine sehr gute Schatzung der fur 

2 5 die realen Daten zutreffenden Dichtef unktion der 

Nullhypothese dar, aus der der gesuchte kritische Testwert 
abgelesen werden kann. 

Als Beispiel ist in Abbildung 2 die mittels der hier 
beschriebenen Methode berechnete Dichtef unktion der 
30 Nullhypothese fiir die Test-Modif ikation 1 angegeben. 

Die Anwendung der erf indungsgemaBen Losung hat folgenden 
Vorteil : 


Wahrend die bekannten Losungen durch Beschrankung auf die 
Ranks der Phasenwinkel (q- sample uniform scores Test) bzw. 
die Ranks der Phasenwinkel und der spektralen Amplituden 

(Modif ikation 1) die im Spektrum enthaltene Information nur 
teilweise nutzen, werden durch die erf indungsgemaiien 
Modif ikationen (Modif ikation 2-4) mehr Information (Mod. 2 
und 3) bzw. die gesamte im Spektrum enthaltene Information 

(Mod. 4) genutzt. Daraus resultiert eine hohere Testpower der 
vorgeschlagenen Modif ikationen . Eine hohere Testpower fuhrt 
dazu, daJi bei einem Horscreening mit vorgegebenen Reizpegel 
die Antwort schneller detektiert wird, dadurch ist der 
Zeitaufwand fur das Screening geringer. Bei einer objektiven 
Horschwellenbestimmung ist infolge der hoheren Testpower eine 
genauere objektive Schwellenbestimmung moglich^ da der 
Antwortnachweis naher an der Horschwelle des Patienten 
gelingt . 

Der Verzicht auf die VerteilungsunabhSngigkeit bei den 
vorgeschlagenen Modif ikationen ist kein Nachteil, da die 
zutreffenden Verteilungen der Nullhypothese entsprechend dem 
vorgeschlagenen Verfahren an Hand der Daten ermittelt werden. 
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Patent anspriiche 

1, Statistische Testverf ahren zum objektaven Nachweis von 
Auditory Steady-State Responses (ASSR) im 

5 Frequenzbereich unter verwendung q- sample uniform scores 

Test/ wobei nur die Phasenwinkeln verwendet werden. 

2, Statistische testverf ahren nach anspruch 1, wobei den 
iiber eine Four iertransf ormation berechneten 

10 Phasenwinkeln verwendet werden . 

3, Statistische Testverf ahren zum objektiven Nachweis von 
Auditory Steady-State Responses (ASSR) im 

Frequenzbereich unter verwendung q-sample uniform scores 
15 Test, wobei spektrale Amplituden und Phasenwinkel 

verwendet/ die Phasenwinkel aber ungerankt, wahrend fur 
die spektralen Amplituden weiterhin die Ranks in den 
Test eingehen. 

20 4. Statistische Testverf ahren zum objektiven Nachweis von 
Auditory Steady-State Responses (ASSR) im 

Frequenzbereich unter verwendung q-sample uniform scores 
Test wobei die mittels der Fouriertransf omation 
berechneten Werte der Phasenwinkel als auch der 
25 spektralen Amplituden direkt (ungerankt) verwendet 

werden . 

5. Statistische Testverf ahren zum objektiven Nachweis von 
Auditory Steady-State Responses (ASSR) im 
30 Frequenzbereich unter verwendung q-sample uniform scores 

Test nach einer den oben genanten anspruche wobei den 
testverf ahren teil eine computerprogram ist, welche auf 
einer lagermedium, sowie einer diskette, einer CD-ROM, 
einer Hard-Disk oder ahnliches gelagert ist. 
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6. Testgerat fur die ausfuhrung einer Statistische 

Testverf ahren zum objektiven Nachweis von Auditory 
Steady-State Responses (ASSR) im Frequenzbereich unter 
5 verwendung q- sample uniform scores Test, mit einer 

computerprogram umfassend eine oder mehrere von den 
obigen genanten f unktionalitaten. 


13 


Zusammenf assung 

Die Erfindung betrift einer statistische Testverf ahren zum 
objektiven Nachweis von Auditory Steady- State Responses 
5 (ASSR) im Frequenzbereich unter verwendung q- sample uniform 
scores Test, wobei nur die Phasenwinkeln verwendet werden. In 

einer Ausf uhrungsbeispiel werden den iiber eine 

Fouriertransf ormation berechneten Phasenwinkeln verwendet. In 
einer anderen Ausf uhrungsbeispiel werden spektrale Amplituden 

10 und Phasenwinkel verwendet, die Phasenwinkel aber ungerankt, 
wahrend fur die spektralen Amplituden weiterhin die Ranks in 
den Test eingehen. In einer noch anderen Ausf uhrungsbeispiel 
werden die mittels der Fouriertransf omation berechneten Werte 
der Phasenwinkel als auch der spektralen Amplituden direkt 

15 (ungerankt) verwendet werden. DieErfindung betrift noch ein 
Testgerat fur die ausfiihrung die Statistische Testverf ahren . 
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